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ABSTRACT 


The flow fields from a turbulent channel simulation are used to compute the budgets for 
the turbulence kinetic energy (k) and its dissipation rate (e). Data from boundary layer 
simulations are used to analyze the dependence of the eddy-viscosity damping-function on 
the Reynolds number and the distance from the wall. The computed budgets are used 
to test existing near-wall turbulence models of the k-e type. We find that the turbulent 
transport models should be modified in the vicinity of the wall. We also find that existing 
models for the different terms in the e-budget are adequate in the region away from the 
wall, but need modification near the wall. The channel flow is computed using a k-e model 
with an eddy- viscosity damping function from the data and no damping functions in the 
e-equation. These computations show that the fc-profile can be adequately predicted, but 
to correctly predict the e-profile, damping functions in the e-equation are needed. 

1. INTRODUCTION 

Of the models used to predict turbulent flows, a popular model is the two-equation 
k-e model. This model gained popularity with the advent of large computers and has been 
used extensively for engineering flows even though it fails to correctly predict a number 
of flows (e.g., the Stanford AFOSR-HTTM conference proceedings.) The most commonly 
used model was developed for high-Reynolds-number flows and is used in wall-bounded 
flows in conjunction with wall functions to patch the core region of the flow to the wall 
region (for a review see, Rodi, 1980). 

The usefulness of a k-e turbulence model that would be valid all the way to the 
wall has been recognized early in the use of the model. Widely used models that use wall 
corrections to the standard k-e model are the models of Jones and Launder (1972, hereafter 
J&L) and Chien (1982, hereafter CH). Several other models with wall corrections have been 
developed and tested for boundary layer flows, each with proposals for damping functions 
and modifications to the high-Reynolds-number k-e model. Bernard (1986) evaluated 
four models by computing the channel flow with these models and found that all models 
underpredict the peak in the turbulent kinetic energy characteristic of near-wall flows. 
He attributed the failure of the models to poor modeling of the pressure-diffusion term. 
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His conclusion is based on comparisons with extrapolations from experimental data in the 
y + < 10 region. We shall see that the pressure diffusion term in the turbulent kinetic- 
energy budget is small compared to the other terms at all y + . In a systematic study, Patel, 
Rodi, and Scheuerer (1985) evaluated eight models for boundary layer flow predictions. 
They found that three models which are based on the k-c model give similar results; the 
model of Lam and Bremliorst (1981) with elaborate damping functions was not superior to 
the CH model or a slightly modified J&L model. They also state that further refinement 
to the existing models is needed for any of the models to be used with confidence. The 
suspected weakness of the model is believed to be the closure models used in the e-equation. 
The current trend in modeling the 6-equation is to introduce ad-hoc damping functions in 
the models of the different terms. An exact assessment of the models has not been possible 
because it is difficult to accurately measure the dissipation rate of the turbulent kinetic 
energy near the wall. 

With the advent of both advanced numerical methods and large-scale computers, 
full simulations of turbulent flows at low Reynolds numbers (Re) are now possible. Of 
interest to the near- wall k-e modeling are the recent simulations of Moser and Moin (1987) 
of a curved channel flow, Spalart (1986a, b) of a flow over a flat plate, and Kim, Moin, 
and Moser (1987, hereafter KMM) of a channel flow. In the calculation of KMM, the 
turbulent flow field for a channel flow at Re — 3300 (based on half the channel width 
and centerline velocity) has been computed and the results compared with experimental 
data. The physical realism of these computed turbulent flowfields has been validated by 
comparing with measured turbulence statistics as well as turbulence structural information 
obtained from laboratory experiments. 

The data base of KMM is used to compute the budget of the turbulence kinetic 
energy (lb), and the budget of the dissipation rate of k (c). These budgets are used to 
test closures for these equations. The data of Spalart are used to analyze the dependence 
of the eddy-viscosity damping function on y + and the Reynolds number. Finally, results 
using an eddy-viscosity damping function with no damping functions in the e-equation 
show that the Jb-profile can be adequately predicted. However, to predict the e-profile, 
damping functions in the e-equation are needed. 

2. THE AVERAGED EQUATIONS 
2.1 The mean momentum equation. 

The averaged Navier-Stokes equations for incompressible flows nondimensionalized 
with the wall variables u T ( = \/t w / p, the friction velocity) and u (the kinematic viscosity), 
are written as follows: 


Vi,t + {ViVjii + («>' L = -P,i/p+u itii (1) 

Ui,i = 0 

where U{ and u\ represent mean and fluctuating velocities respectively, and P is the mean 
pressure. In most simple phenomenological models, a Boussinesq approximation is used to 
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close Eq.(l). In this case, the terms u\u'- are approximated as 


u^Uj — 2u r pSij 


( 2 ) 


where k = u\u’J 2 is the turbulent kinetic energy, - (t/jj + Uj,i )/ 2 is the mean strain- 
rate tensor, vt is an eddy viscosity, and f> t j is the Kronecker delta function. 

The k-e modelers approximate the eddy viscosity vj as follows (see, for example, Patel 

et «/., 1985): 

i'T = Cufp— ( 3 ) 

e 

where C ^ is a model constant, is a damping function, and l is a modified dissipation 
rate of k. Note that e = u' , + D where D is a function that has the same dimension 

l 'J l 'J 

as t = u\ ju\ it is chosen so that t vanishes at the wall and e — t at high-Reynolds 
numbers. Different near-wall fc-e models use different expressions for D. Note that in this 
approach, the problem was reduced from deriving equations describing the evolution of 
u\u' } to deriving equations for k and e. 

2.2 The ^-Equation 

For an incompressible flow the exact equation governing the transport of k is 

k, t + {Ujk),j = P k + T k + n fc + D k - e (4) 

where the different terms on the right-hand side are given as (rate of ...), 


Pk = -U^u'jSij 

Tk = 

=--(u' l p‘), t 

P 

Dk —^,jj 


e 




Production 
Turbulence transport 
Pressure diffusion 
Viscous diffusion 
Dissipation 


Figure 1 shows the various terms computed from the channel data of KMM as a func- 
tion of the wall variable y + . The salient feature of this plot is that away from the wall, the 
production rate is almost balanced by the dissipation rate. Close to the wall, the produc- 
tion rate and dissipation rate are still the dominant terms, but the turbulent transport rate 
and viscous diffusion rate are no longer negligible. Only the pressure diffusion rate remains 
negligible compared to the other terms. The turbulent diffusion rate has a positive peak 
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a t y + — 6 and a negative peak at — 15. The viscous diffusion rate and the dissipation 
rate are related by the identity 


£ = ~ u i u i,jj = e " (5) 

At the wall, the left-hand side of Eq.(5) vanishes and the diffusion rate exactly balances 
the dissipation rate. The expansion of the left-hand side of Eq. (5) in Taylor series will 
yield (' = 0(y+) in the vicinity of the wall and according to the data, it has a positive 
slope at the wall. Away from the wall, the second derivative of k is small compared to (or 
of the opposite sign as) e which implies that e' is positive throughout the channel. In fact 
for homogeneous flows, e' = e. 

Comparison of the budget data in Fig. 1 with the estimates using Laufer’s data (see 
Townsend, 1980, p. 145, Fig. 5.5) shows that both the turbulence transport and pressure 
diffusion terms are overestimated. The viscous diffusion at the wall is underestimated in 
Laufer’s data which yields a lower dissipation rate at the wall. Our data are consistent 
with the budget data of Moser and Moin (1987) for a flow in a curved channel and of 
Spalart (1986a, b) for flows over a flat plate; close to the wall, all simulation data show 
that the pressure diffusion term remains small compared to the other terms in the budget. 

2.3 The f-Equation 

The equation describing the evolution of e is 


e ,t + (Vje)j = Pi + Pi + Pi + Pl + T ( + n ( + D ( -r 


(6) 


We can identify the different terms on the right-hand side as (rate of ...), 


p; =-2 




P 4 — — f rU 1 UV 
f i, k i,m k ,m 


t { = ~K< m < m ),* 


n e =--(p[mK,m),k 

D( =e,fcfc 


T = 2 u[ 




Production by mean velocity 
gradient 

Mixed production 
Gradient production 
Turbulent production 
Turbulent transport 

Pressure transport 
Viscous diffusion 
Dissipation 


Tennekes and Lumley (1972) analyzed the vorticity fluctuation budget which is related 
to the above budget for homogeneous flows. They inferred from an order-of-magnitude 
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analysis that in the high-Reynolds-number regime, the turbulent-production rate (P*) 
and dissipation rate (T) dominate the balance equation. However, the difference of these 
terms yields a term of the order of the other terms. The various terms in the balance 
equation for e are shown in Fig. 2. The present results indicate that P 4 and T are the 
largest terms in the core region of the channel in agreement with the analysis of Tennekes 
and Lumley. Near the wall, these terms are still large but are not larger than the other 
terms. Close to the wall (y + < 8), the production rate P] becomes of the same order as 
P 4 . In the range 6 < y + <15, the mixed-production rate (P 2 ) is of the same order as P 4 . 

3. TURBULENCE MODELS 


In the previous section we introduced the eddy-viscosity model to represent the 
Reynolds stresses and in turn modeled the eddy viscosity in terms of a damping func- 
tion, the turbulence kinetic energy (fc), and the dissipation rate of k («). In turbulence 
modeling using k and e, all the terms in the balance equations for k and e that involve 
correlations other than k and e have to be modeled in terms of Ui, k, and e. In this section, 
we use the simulation data to show that the eddy- viscosity damping function should be a 
function of both the Reynolds number and y + . We also use the simulation data to test 
closure models for the k and c equations. 


3.1 Damping the eddy viscosity 

In flows where the relevant gradients are in one direction only (e.g., a fully developed 
channel or homogeneous shear), the Boussinesq approximation and the definition of the 
eddy viscosity will yield the following: 


u' t u 2 

~u7 


k 2 


( 7 ) 


It can be shown that u\u 2 = 0(y + ), and k = 0(y + ) as y + *—► 0. The choice of the 
near-wall behavior of the damping function, / M , will depend on the near-wall behavior of 
l. Using D = -kjj will yield i = (' = 0(y + ) in the vicinity of the wall. 


The equation for the production rate of the turbulence kinetic energy (Pk = — ain't/, 2 ) 
can be combined with the definition of the eddy viscosity and the Boussinesq approximation 
to yield 


c,u = 


■(¥) 


Pk 


(8) 


The damping function, / M , should be constructed by examining the behavior of the 
individual terms -u\u 2 /k and P*/e' in wall-bounded flows. Figure 3 shows the distribution 
of -u\u' 2 /k as function of y + from the bo undar y layer data of Spalart (1986b) for three 
Reynolds numbers. The data show that -u\u' 2 /k is a function of the Reynolds number, 
and for large y + , the term seems to asymptote to a value close to 0.3, 

-u\u 2 /k = /«v(y + ,Pe) (9) 
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The variation of Pfc/e' as a function of y + was also computed from the data of Spalart 
(1986) and KMM (see, Fig. 4). In the vicinity of the wall (y + < 40), the data show a 
reasonable collapse and therefore exhibit dependence only on y + , 

P*/<’ = f r/ ,{y + ) (io) 

The function f p / t peaks at around y + = 19 and f p / f = (?(y +2 ) in the vicinity of the wall. 

3.2 ^-Balance 

In the ^-equation, the viscous diffusion term need not be modeled, but the pressure- 
diffusion term and turbulent-transport term are usually added and modeled as one term, 

T k + 11* = {v T kj) %j (11) 

Figure 5 shows the distributio n of t he turbulence transport term compared to the eddy- 
viscosity model using ut = ~u\ u^/Uj from the data. In the vicinity of the wall, the 
model has a different slope than the data would indicate; this can also be shown from a 
Taylor-series expansion in the vicinity of the wall. 

The production term is modeled by substituting the model for — uju' in the production 
rate expression, 

-u'iu'jUij - u T 2 SijSij ( 12 ) 

In this case the data are matched exactly. 

3.3 e-Balance 


The e-equation is closed by modeling the te rms in Eq. 6 . To r epresent the first 
two production terms we need an expression for u' i k u '- k and u' k i u' k ■ in terms of it, e 
and other mean quantities. Note that these two terms have the same trace; they are 
related for homogeneous flows through the vorticity fluctuation. If we assume that Rotta’s 
approximation is valid for both terms, 


u'u' 

l J 


+ <.X,,) = -ft -r 2 ' 


(13) 


and substitute the Boussinesq approximation for the Reynolds stresses we obtain: 


Pi + P? = C i ^v T 2S ij S ij (14) 

This expression is the same as the commonly used model for the production of c; Figure 6 
shows that model compared to the exact expression, using the constant recommended by 
CH, Ci = 1.35. The model yields a lower peak than the data would indicate. For near-wall 
models, the common approach (see Patel et ai., 1985) is to introduce a damping function, 


fli 


1 + 2 


_ P[±P[ 

C\C p k2S t} Sij 


(15) 
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The term P 2 is negligible compared to the other terms in the channel flow, so no 
explicit expression is used to model it. We expect the turbulent production term (P 4 ) 
to be non-negligible even in isotropic flows; an appropriate model for this term will be 
a function of k and e and not a function of the mean velocity. Dimensional analysis 
yields P 4 oc e 2 /k. The same arguments are used in modeling the dissipation rate of e as 
oc t 1 jk. The two terms are then combined and modeled propotional to f 2 /k. However, 
near the wall t 2 jk i— > oc and the term should be modified for near-wall effects. This is 

achieved by using the modified dissipation (e" - t - 2^(fc 1/2 ) 2 ) ) of J&L and setting 

P 4 — T oc ee" jk. In this case e" jk and e are bounded as y + *— ► 0. Hanjalic and Launder 
(1976) proposed a similar model and from experimental data of grid turbulence, they 
inferred that the proportionality factor should be a function of the turbulence Reynolds 
number {Re t = k 2 / c); therefore, they added a damping function of Rc t . We found (by 
examining the channel data) that this damping function is an unnecessary complication 
to the model. The model for the combined terms is given as 


P 4 -T = -C 2 ^ (16) 

Figure 7 shows the comparison of this model with the data and with the model of Hanjalic 
and Launder (6' 2 = 1.8 in their model). The difference between the models is smaller than 
the difference between the data and the models. The models adequately compare with the 
data in the y + > 11 range, but underpredict the data close to the wall. We have seen in 
Fig. 6 that the production rate is underpredicted in the < 15 region and as we will see 
from the computational results, the underprediction of both the dissipation rate and the 
production rate are the major source of discrepancy between the predictions and the data. 
The common approach in near-wall models is to introduce in conjunction with a damping 
function for the production rate of e, a damping function for the dissiaption rate, 


T - P 4 


(17) 


The remaining terms in the balance equation of e are transport rate terms, which are 
grouped together and modeled using an eddy- viscosity-diffusion model, 

T t + n ( = (— f,j) (is) 

where a — 1.3 (J&L and CH). Figure 8 compares the model with the terms representing 
the left-hand side of Eq. 18. As in the case for the fc-equation, the comparison is not good 
in the vicinity of the wall, in fact, it can be shown from Taylor series expansion that this 
model does not have the proper asymptotic behavior as y + >-+ 0. 

4. RESULTS USING k-e MODELS 

In our model-testing in the previous section, we assumed that we know the distribution 
of U , k, and e; we then tested the model expression. If the models agreed perfectly with 
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the data for all the terms and if the resulting set of equations admit a unique solution, 
we can expect the computational results to yield accurate predictions. In general, the 
agreement is not perfect. For example, the production rate and dissipation rate of e are 
not well modeled in the < 15 region; therefore we have no assurance that the model 
will yield accurate predictions. In this section, we will compute the channel flow using 
the k-t model described in the previous section using the eddy-viscosity-damping function 
from the data and no damping functions in the c-equation. 

In our prediction of the channel flow we begin with the averaged Navier-Stokes equa- 
tions and simplify them for the case of a fully developed channel. The boundary conditions 
used are U — 0, k = 0, e = k >2 2 at the wall, and f/ 2 = k <2 — (,2 = 0 at the centerline 
(y + = Re r — 180). The equations were solved on a nonuniform mesh using the fourth- 
order Runga-Kutta method to discretize the equations in time and central differencing for 
space discretization. The same numerical scheme was used to carry out three calculations 
using the k-t models of J&L, CH, and the model described in section 3. 

The mean velocity profiles as predicted by all models give acceptable results but using 
the exact damping function gives a better overall agreement with the data. Figure 9 shows 
the Ar-profiles as predicted by the near- wall models. The model of J&L does not reproduce 
the sharp peak in the kinetic energy expected near the wall. The model recommended 
by CH predicts a peak at a slightly shifted location, but the level of k is overpredicted in 
the t/ + > 20 range. When the damping function from the data is used, the peak level is 
reproduced at approximately the same y + . The level of k is also better predicted in the 
rest of the channel. These predictions indicate that the eddy-viscosity damping function 
plays a key role in the prediction of the fc-profile. 

Figure 10 shows the profiles of t from the three computations. We find that the models 
of CH and J&L predict the same total-dissipation-rate e. Near the wall, all models fail to 
predict t correctly. The models predict a peak in t at y + ss 10 whereas the data shows 
that t peaks at the wall. However, it is interesting to note that the peak predicted by 
the models is at the same location where the data show a local peak. This comparison 
indicates that the damping functions in the e-equation will play an important role in our 
prediction of the correct e-profile. 

5. CONCLUSIONS AND DISCUSSION 

In this paper we used turbulence data from a full simulation to compute the terms in 
the budgets of the turbulent kinetic energy and of the dissipation rate. 

These budgets were used to test closure models for the k and c equations. We find that 
the transport models need to be improved near the wall. The closures to the e-equation 
are good away from the wall but poor in the vicinity of the wall. We have split the eddy- 
viscosity damping function into two functions; one representing the ratio of the production 
over a modi fied e which we find to be function of j/ + only, the other representing the 
parameter -u\u' 2 /k which we find to be a function of t/ + and the Reynolds number. 

Computations with the k-t model indicate that the eddy- viscosity damping function 
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plays a key role in the prediction of the fc-profile. To improve the prediction of the e-profile, 
damping functions in the e-equation need to be introduced. In fact, using the damping 
functions / ( \ +2 and f ( 2 from the data gives a better prediction for the e profile (see Fig. 
11). However, the improvement is not complete. By neglecting P 3 , we effectively have 
lumped the term with the model of P} + P? and the damping function for the production 
term should be 


ri + 2+3 _ 
J ( 1 — 


P] + P; + P't 

C,C M fc2SySij 


(19) 


Use of the above damping function yields a better c-profile (see Fig. 11). The remainder 
of the disagreement must come from the transport terms. 


ACKNOWLEDGEMENTS 


We are grateful to Dr. P. Spalart for generously providing us with the data used to 
develop the eddy- viscosity damping function. We are also grateful to Drs. M. J. Lee, 
T. H. Shih, and Prof. W. C. Reynolds for numerous discussions we had during the course 
of this work 


REFERENCES 

1980-81 AFOSR-HTTM-STANFORD Conference on Complex Turbulent Flows, 
Vols. II and III, S. J. Kline, B. J. Cantwell, and G. M. Lilley eds., Thermosciences Division, 
Mechanical Engineering Dept., Stanford Univ., Stanford, Calif., 1982. 

Bernard, P. S., “Limitations of the Near- Wall k-e Turbulence Model,” AIAA J., Vol. 
24, No. 4, pp. 619-622, 1986. 

Chien, K. Y., “Predictions of Channel and Boundary-Layer Flows with a Low 
Reynolds Number Turbulence Model,” AIAA J., Vol. 20, No. 1, pp. 33-38, 1982. 

Hanjalic, K., and Launder, B. E., “Contribution Towards a Reynolds-Stress Closure 
for Low-Reynolds-Number Turbulence,” J. Fluid Mech., Vol. 74, pp. 593-610, 1976. 

Jones, W. P., and Launder, B. E., “The Prediction of Laminarization with a Two- 
Equation Model of Turbulence,” Int. J. Heat Mass Transfer, Vol. 15, pp. 301-314, 1972. 

Kim, J., Moin, P., and Moser, R., “Turbulence Statistics in Fully Developed Channel 
Flow at Low Reynolds Number,” J.F.M., Vol. 177, pp. 133—166, 1987. 

Lam, C. K. G., Bremhorst, K., “A Modified Form of the k-e Model for Predicting 
Wall Turbulence,” Trans. ASME, J. Fluids Eng., Vol. 103, pp. 456-460, 1981. 

Moser, R. D., and Moin, P., “The effects of Curvature in Wall-Bounded Turbulent 
Flows,” J.F.M., Vol. 173, pp.479-510, 1987. 

Patel, V. C., Rodi, W., and Scheuerer, G., “Turbulence Models for Near-Wall and 
Low Reynolds Number Flow: A Review,” AIAA J. Vol. 33, No. 9, pp. 1308-1318, 1985. 


9 



Rodi, W., “Turbulence Models and their Application in Hydraulics,” International 
Association of Hydraulic Research, Monograph, Delft, the Netherlands, 1980. 

Spalart, P. R. Numerical Study of Sink-Flow Boundary Layers. J. Fluid Mech., Vol. 
172, pp. 307-328, 1986a. 

Spalart, P. R. Direct Simulation of a Turbulent Boundary Layer up to R # — 1400,” 
NASA TM 89407. NASA-Ames Research Center, Moffett Field, CA, 1986b. 

Tennekes, H., and Lumley, J. L., A First Course in Turbulence, The MIT Press, 
Cambridge, Mass. 1972. 

Townsend, A. A. 1976 The structure of turbulent shear flow. Cambridge Uni- 
versity Press, Cambridge, England. 


10 




0.03 


0.02 

Gain 

0.01 

0.00 

- 0.01 

Loss 

- 0.02 



0.0 20.0 40.0 60.0 



Figure 2. Terms in the budget of the dissipation rate of the turbulence kinetic 
energy, t, in wall coordinates. P] = production by mean velocity gradient; P] = mixed 
production; P ( 3 = gradient production; P* = turbulent production; T ( = turbulent trans- 
port; D ( = viscous diffusion; T = dissipation rate; II ( = pressure transport. 
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Figure 3. Distribution of — u^u^jk in wall-bounded flows. Boundary-layer data is 

from Spalart (1986b). 
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Figure 4. Distribution of P*/e' in wall bounded flows. Channel data is from KMM 
and boundary-layer data is from Spalart (1986b). 
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Figure 5. Transport terms Tk + life across the channel, o o o o term computed 
from the channel data; model, Eq. 11. 
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Figure 9. Turbulence kinetic energy distribution for Re T = 180. □ □ □ □ term 

puted from the channel data; model using / M from the data, Chien’s 

lei. Jones &; Launder’s model. 



Figure 10. Dissipation rate (e) distribution for Re r = 180. o □ □ □ term com 

puted from the channel data; model using from the data, Chien’s model 

Jones & Launder’s model. 



Figure 11. Dissipation rate (e) distribution for Re r = 180. □ □ □ □ term com 

puted from the channel data; model using from the data with f\ j +2 = f ( 2 = 1. 

model with / M , /*, +2 , and / t2 from the data. model with / M , f] i +2+3 , and f (2 

from the data. 
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